Population Ecology

Demographic models are crucial for uncovering the mechanisms underlying changing populations trajectories and structure and to identify the drivers of such changes. Common pochard populations of North east/North west Europe have experienced a sharp decline accompanied by an increase in male proportions among adults since the mid-1990. We used a two-sex, two-stage deterministic matrix population model, to perform a prospective perturbation analysis and to explore, through simulations, some plausible causes underlying the observed decline and change in sex structure. We show that Common pochard populations are more sensitive to changing survival than to changes in productivitys components (clutch size, nest survival...). However, due to an environmental variance much higher than that of survival, components of productivity, especially nest survival, would be the main drivers of Common pochard populations growth rate, a finding supported by empirical data. More importantly, we show that although sex-specific changes in survival at any stage of the life cycle are potent drivers of both population growth rate and changing sex ratio, there is no need to resort to them for explaining the increasing proportions of males such as observed in Common pochard. Because adult males display higher survival than adult females (0.74 against 0.64 on average), any factor affecting recruitment (nest or first year survival) increases the weight of adults into the populations and hence the proportion of males. Thus, in species displaying sex-biased mortality, such as many ducks, decreasing recruitment can underly declining population size and changes in sex structure at the same time, emphasising the importance of accounting for males in monitoring schemes and demographic models.

Arctic and Boreal raptor communities will continue to be affected by borealization and other climate change related processes, providing a challenge for ecologists predicting future sates. However, by using community assembly theory and species traits, future communities may be predictable. In this study, we analyzed variation in reproduction traits as a consequence of diet specialization for 29 raptors, 2 skuas and 3 corvids. We assessed and implemented foraging traits for specialists and generalists into predator-prey models from which successful invasion conditions were derived. Specialist raptors produced larger clutch sizes, had a higher proportion of fledged per clutch and also expressed more variation compared to generalist raptors. These results suggest a relationship between diet specialization and reproductive traits which was also observed within phylogenetic orders. Specialist owls (Strigiformes) produced higher clutch sizes with a larger clutch range compared to generalist owls. The same pattern was observed for falcons (Falconiformes). No clear difference in reproduction was observed for specialist and generalist hawks, kites and eagles (Accipitriformes). Corvids expressed clutch sizes similar to that of specialist raptors while having the lowest proportion of fledged per clutch. Differences in foraging traits between specialists and generalists could be distinguished using functional response curves. A predator-prey model parameterized with foraging trait data showed that a generalist can coexist with a resident specialist if it has access to prey unavailable to the resident specialist. Otherwise, the native specialist outcompetes the invading generalist due to foraging efficiency. The combined empirical and theoretical findings in this study show how diet specialization affects both reproduction and the potential invasion success of raptors.

Age-dependent survival is central to understanding population dynamics and life-history evolution. We analysed carcass weight and age-at-harvest data from 6022 red foxes (Vulpes vulpes) collected across Sweden between 1967 and 1971 to evaluate latitudinal effects on body mass and age-dependent survival. Carcass weights decreased from south to north in both adults and sub-adults, contrary to Bergmann's rule, with southern foxes weighing approximately 1.27 times more than northern foxes. The latitudinal weight gradient exceeded the sex difference in both age classes, and no sex x region interaction was detected. The decrease in weight with latitude is consistent with reduced prey availability and harsher winter conditions in the north, which limit growth and body size during development. Using a Bayesian age-at-harvest model with region-specific population growth rates (lambda), we estimated age-dependent survival probabilities for four latitudinal regions and both sexes. Despite the strong latitudinal gradient in weight, survival did not show a corresponding pattern - regional differences were uncertain, with all credible intervals spanning zero. Regional population growth rates were consistent with slight decline in the north and near-stability in the south-central region, which suggests that body condition and population dynamics are coupled at the regional scale despite no survival gradient. The decoupling of body condition and survival across regions suggests that mortality patterns are similar across the latitudinal gradient. We discuss these patterns in terms of latitudinal productivity gradients, prey availability, and life-history trade-offs in a widely distributed carnivore.

Humans are efficient and deadly predators, yet they may also interact with wildlife in non-lethal ways. This study explores how interactions with lethal and non-lethal human "superpredators" alter predator-prey dynamics using an agent-based modelling approach. Our model incorporates both the consumptive (lethal) and non-consumptive (behavioural) effects of humans, as well as of predators on prey. We explored how the replacement of apex predators by humans affects mesopredator-prey dynamics, with particular emphasis on trophic targeting and differences between lethal and non-lethal interactions. We found that human superpredators have a greater effect on model outcomes than apex predators. When superpredators consume mesopredators alone or with prey, the probability of mesopredator-prey coexistence increases to a greater extent than when apex predators consume mesopredators. In contrast, superpredators consuming only prey slightly increases overall extinction risks and reduces coexistence. Non-lethal superpredators, despite eliciting anti-predator responses in mesopredators and prey, had a negligible effect on population dynamics. Our findings demonstrate that human superpredators may functionally replace apex predators when they are lethal. However, non-lethal interactions with humans may not be as ecologically significant as lethal interactions, even when humans induce anti-predator responses.

Environmental variability can strongly alter coexistence among competing species and their extinction risk, particularly when population dynamics are shaped by behavioral interactions, such as fear. In this work, we develop a novel stochastic differential equation competition model that incorporates both non-consumptive fear effects and environmental variability to investigate how behavioral interactions influence species coexistence under random fluctuations. Our result reveals that environmental stochasticity can drive species to extinction even when the corresponding deterministic system admits coexistence. In particular, under an explicit stability condition on the fear and competition parameters and sufficiently strong averaged noise intensities, we prove that both competing species become extinct exponentially almost surely. Conversely, we derive a stochastic persistence criterion in terms of fear, competition, and noise-induced suppression parameters for the fearful species. We further demonstrate that environmental noise may reverse classical competition-exclusion outcomes, leading to qualitatively different long-term dynamics from those predicted deterministically. These results provide rigorous thresholds separating stochastic extinction from persistence and highlight the critical role of environmental variability in fear-mediated competitive ecosystems. From an applied perspective, these results provide insight into how behavioral interactions and environmental variability influence species survival, with potential applications in ecological management and conservation.

Understanding and predicting extinction risk is a central challenge in population biology. Mathematical models incorporating Allee thresholds are commonly used to understand population dynamics and to assess extinction risks. Inaccurate predictions can have serious consequences for conservation management. In this simulation study, we develop a likelihood-based inference and prediction workflow to estimate parameters, including the Allee threshold and population diffusivity parameters, using noisy count data generated using a well-defined discrete model. Although parameters are identifiable according to commonly used criteria, the accuracy of resulting predictions depends strongly on the quantity, quality, collection time and spatial resolution of the data. Our workflow demonstrates that seemingly reliable parameter estimates can lead to inaccurate predictions, highlighting the need for careful consideration of data quality and quantity to guide extinction-risk modelling and prediction. Open source software is provided on GitHub to replicate and extend all results considered.

The mosquito Aedes albopictus is a competent vector for the transmission of several arboviruses and is currently spreading across many continents. Since conventional control methods, like insecticides, often lead to environmental problems and the emergence of resistance, scientists developed alternative mosquito control strategies. One of the most used is the Sterile Insect Technique (SIT), which involves the mass release of males sterilized through irradiation. The Toxic Male Technique (TMT) is instead based on the release of genetically modified males expressing toxic proteins that kill females when they mate. Control strategies are often intended as methods to eradicate mosquito populations, yet a less ambitious and more cost-effective task is to reduce them such that the probability of transmission of viruses to humans becomes negligible. To compare the efficacy of these control strategies, we develop a mathematical model with two communicating compartments: a mosquito population and epidemiological model coupled with a human epidemiological model. As a proof-of-concept, we test the model using meteorological and entomological data for the Emilia-Romagna region. Our results indicate that the TMT strategy is more effective in lowering the probability of transmission and provides indication for the deployment of control strategies.

In order to feed a growing global population without silencing nature, conceiving agricultural management strategies reconciling yield and conservation goals is key. Using numerical simulations of a metacommunity model, we explore the possibilities for compromise offered by spatial management strategies of farmed areas. Each strategy is characterized by its farming intensity, the proportion of farmed lands and their spatial aggregation. We show that achieving equitable yield-biodiversity compromise is difficult. While conciliatory strategies offering top yield and biodiversity are typically not possible, accepting slightly lower yields (ie, "Pretty Good Yield strategies") allows to recover substantial biodiversity. Such reconciliation possibilities are limited for species with small dispersal. Yield increases mainly through farmland expansion, whereas farming intensity strongly influences biodiversity, increasing it at low intensity before decreasing with further intensification. Finally, we demonstrate that reconciliation is easier if agricultural production relies on biodiversity through ecosystem services.

Rapid urbanization has intensified ecological degradation and threatened ecosystem sustainability in the Yangtze River Basin. This study proposes a Function-Structure-Dynamic (FSD) framework integrating ecological protection importance assessment, ecological security pattern (ESP) construction, and future land-use simulation to identify ecological restoration priorities in the Anhui section of the Yangtze River Basin. Ecologically important areas were identified by coupling ecosystem services and ecological vulnerability. MSPA, landscape connectivity analysis, resistance surface modeling, and circuit theory were applied to construct the ESP, while the PLUS model simulated land-use change and ecological risks in 2032. Results identified 37 ecological sources covering 6,901.39 km{superscript 2} and 84 ecological corridors with a total length of 916.80 km. Ecological pinch points, barrier areas, and 17 ecological warning areas were further delineated. Based on these findings, four ecological restoration zones were proposed with differentiated management strategies. The FSD framework effectively integrates ecological functions, spatial structure, and future dynamics, providing scientific support for ecological restoration and sustainable territorial spatial planning in rapidly urbanizing regions.

Downstream migration and sea entry are periods of high mortality for sea trout smolts and migration timing is a critical aspect of survival. We aimed to investigate migration timing of PIT-tagged sea trout smolts across a latitudinal gradient in five freshwater systems across Europe: Norway, Sweden, Denmark, France and Portugal. In two systems, smolt migration was further examined to assess a) differences in sex, body size and condition between autumn and spring migrants; b) the influence of spatial origin within the stream; and c) relationships between individual body size and migration date. Tagged smolts were not detected migrating in the Portuguese watershed. Spring migration timing differed significantly between watersheds in Norway, Sweden, Denmark and France. Generally, there was a trend of earlier migration at lower latitudes. Autumn vs spring migration was examined in Gudso-Denmark. Autumn migrants were larger in both length and mass, with no differences in sex ratios or body condition. Fish originating from an upstream site were more likely to migrate in the autumn compared to the spring and vice versa. Size dependent migration was found in the Swedish system, Haga [a]-Sweden, with larger individuals migrating earlier in the spring than smaller individuals. Outward-migrating smolts were also more likely to originate from a downstream site than an upstream site. Overall, these results show both large-scale geographic and fine-scale individual influences on migration timing. Given that climate change may have large impacts on migration patterns in sea trout, understanding variability in migratory patterns across a latitudinal gradient is an important tool for predicting responses to environmental changes.

Fledging is a critical period of a seabird life cycle. Using satellite telemetry, we compared movements and survival proxies (transmission duration) of chinstrap penguin fledglings tracked in 2017 (n=8) and 2025 (n=17) relative to krill fishing vessel activity. In 2017, fishing vessels operated intensively near colonies during summer, resulting in early, frequent encounters (median 1.3 days post-fledging) and short transmission durations (median 9.2 days). In 2025, reduced fishing delayed encounters (median 10.0 days) and tripled tracking duration (median 24.0 days). Hidden Markov Models revealed that vessel encounters reduced the probability of transitioning from foraging to transit behavior ({beta} = -0.76), an effect stronger than the productivity ({beta} = -0.11). While 87.5% of 2017 fledglings ceased transmission prematurely within weeks (half of those right after entering areas intensively used by fishing vessels), 65% of 2025 fledglings survived beyond March, with half of those five transmitting until May after dispersing eastward to the South Orkney Islands. These findings suggest that spatiotemporal overlap with krill fisheries during the critical post-fledging window affected foraging behavior and was associated with shorter transmission durations. Our results support further research of post-fledging penguin ecology to better understand the potential impact of fishery, and, following the precautionary principle, support fishing seasonal protection of important areas during critical periods of krill predators life cycle.

O_LIWhen eradication is unfeasible, invasive predator control should evaluate how removal affects ecological responses by native species. Assessments often use total invasive predator abundance to evaluate prey responses, yet intraspecific variation in diet and space use means that some subgroups cause disproportionate impacts. Identifying these problem individuals, and the spatial scales over which their impacts operate, can enable targeted spatially explicit removal to maximise impact reduction. However, despite individual-level information is often already collected during trapping operations it is seldom included when analysing predator impacts, potentially biasing the conservation outcomes expected under blanket removal. C_LIO_LIWe use a novel framework and two decades of invasive predator control data to estimate how individual variation in residency status influences the distance-dependent impacts of invasive American mink Neogale vison on water vole Arvicola amphibius occupancy across two prey surveys. We also develop a sub-model to predict mink residency status for individuals with missing age data. C_LIO_LIThe probability of capturing adult mink decreased with elevation and years of control, indicating that long-term control altered the resident population and demographic composition of mink around water vole sites. C_LIO_LIDistance-dependent negative impacts of mink varied by residency status, becoming negligible at approximately 20 km from water vole sites for resident mink and 2 km for transient. The spatial scale of mink impacts was largest during the first vole survey when resident mink were more abundant, and declined rapidly for the second survey, when mink were less abundant and spatially clustered. Our results suggest that water voles have benefited mostly from reducing resident mink rather than the total population, especially in early control phases. C_LIO_LIManagers can use our framework to develop spatially explicit and impact-based strategies, not restricted to invasive species control, to construct empirically informed management buffers around populations of conservation concern. Long-term efforts will change the landscape and invasive predator contexts, and thus we recommend iteratively updating and re-evaluating management outcome evaluations. We argue that incorporating individual heterogeneity improves our understanding of ecological mechanisms influencing management success but that the suitability of targeted strategies should be evaluated for target socioecological contexts. C_LI

Aims: Species' distributions are determined by the interplay between ecological niche and dispersal ability, constrained by biogeographical barriers. Bald-headed uakaris (Cacajao spp.) are highly specialized primates often associated with seasonally flooded forests. In this study, we used ecological niche models to assess changes in habitat suitability and geographic distribution of uakari species under future scenarios. Location: Western Amazonia. Methods: We integrated ecological niche models, current deforestation data, and dispersal ability to estimate habitat suitability under two Shared Socioeconomic Pathway (SSP) scenarios: intermediate (SSP2-4.5) and very high (SSP5-8.5) greenhouse gas (GHG) emissions. Results: Our models project shifts in suitable conditions for all species. Three of the five species are projected to experience substantial reductions ([≥]62%) in suitable habitat conditions within their current ranges by 2050 under both future scenarios. Across the western Amazonia, up to 219,189 km2 and 211,276 km2 of land are projected to be unsuitable within the uakari ranges under the intermediate and very high emissions scenarios, respectively. This is particularly relevant for C. calvus, C. rubicundus, and C. ucayalii. At the species level, the uakaris may lose between 343 km2 and 84,531 km2 of their ranges in the intermediate scenario and 858 km2 and 76,216 km2 in the very high scenario. Shifts in suitability due to climate change are expected to vary from 6 to 191 km in the intermediate scenario and from 5 to 168 km in the very high scenario. Furthermore, the uakaris may lose between 0.5% and 8% of their current ranges due to deforestation in all scenarios. Main conclusions: Our findings reveal a high sensitivity of the uakaris to climate change impacts. It is projected that all species may experience contractions in the suitable areas and spatial suitability within their ranges by 2050, underscoring climate change as a relevant threat to these taxa.

Fistulobalanus albicostatus Pilsbry, 1916 (Thoracica: Balanidae) is a tropical to temperate species distributed in the NW Pacific. The previously known northernmost record of this species in Japan was from Aomori Prefecture, at the northern end of Honshu Island, Japan. However, field surveys conducted in 2023 and 2026 confirmed the occurrence of F. albicostatus in Hakodate Bay, the southern end of Hokkaido, Japan, across the Tsugaru Strait, thereby extending the northern limit of its known distribution. A line transect survey conducted in May 2026 recorded seven living individuals within an area of 128 m. F. albicostatus was rare on the mid-high shores, accounting for only 2% of all barnacle individuals in a quadrat survey. The basal diameter of the living individuals ranged from 0.76 to 1.23 cm, and all individuals possessed ovaries. Based on characteristics of both morphological and COI gene, the specimens were identified as F. albicostatus, and belonged to the same haplotype of populations that are present in Honshu Island. The establishment of F. albicostatus in Hokkaido suggests an ongoing northward range shift of this warm-water species, with the potential for further expansion under continued ocean warming.

Spatial population genetic and ecological modeling is often necessary to predict outcomes accurately. One example is gene drive, a rapid process involving spread of gene drive alleles through a population, usually to suppress pests or reduce transmission of vector-borne disease. Several existing models have been used to assess gene drive and other spatial processes. However, each of these has limitations, such as high computational cost and limited scalability, difficulty in incorporating environmental factors and complex lifecycles, or potentially simplified spatial structure. To overcome these challenges, we propose a hexagon-based computational framework that is designed to mimic continuous space for rapid genetic wave advances. This allows us to accurately simulate a larger spatial domain with lower computational investment. We implemented this model and compared the wave speeds of different gene drives with those obtained from other models. The results showed good agreement when hexagon width and dispersal were properly calibrated. We then determined optimal circular and linear (along roads) release patterns for a variety of gene drives and Wolbachia bacteria. To demonstrate the application of our framework to a hypothetical scenario, we constructed a model Culex quinquefasciatus mosquitoes on Hainan Island. We then evaluated the outcome of different gene drive release strategies, showing the transgenic insect release level necessary to achieve high gene drive coverage and how this could be further optimized based on mosquito and human distribution. Overall, our hex-based population genetic framework provides a flexible platform for realistic and large-scale models for gene drive and related applications.

Effective conservation depends on demographic metrics that reliably reflect species status, particularly population abundance. For elusive species occurring at low densities, however, such metrics remain difficult to obtain. Spatial capture-recapture (SCR) models are the standardized approach for estimating density in marked populations, but their data requirements, especially the need for multiple spatial recaptures across individuals, often limit applicability in small or data-poor populations. This constraint has resulted in knowledge gaps for some of the most vulnerable species, undermining evidence-based conservation planning and management. Using camera-trap data and SCR-derived density estimates from data-rich populations, we evaluated alternative, less data-demanding metrics and tested the hypothesis: Space to Event (STE), Mean Local Abundance (MLA), and Relative Abundance Index (RAI) exhibit predictable relationships with SCR-derived density; if supported, these metrics can reliably estimate density in populations where SCR models cannot be implemented. We applied this framework to the European wildcat (Felis silvestris), an elusive small felid with highly fragmented populations across Europe, for which density estimates are largely lacking despite growing conservation concern. Across 21 study areas spanning most of the species' range, our results indicate that European wildcats generally occur at lower densities than previously reported. SCR-derived estimates (n=10) averaged 10.32 {+/-} 11.56 inds/100km2, while STE enabled density estimation in five additional data-poor areas (mean 5.52 {+/-} 5.33 inds/100km2). STE showed a strong linear relationship with SCR-derived density (R2=0.98), supporting its use as a viable alternative when SCR is infeasible, although it tended to underestimate compared to SCR, especially at higher densities. In contrast, MLA and RAI showed weaker and non-linear relationships with SCR-derived density (R2=0.65), indicating substantially lower explanatory power and suggesting their estimates are more strongly influenced by confounding processes. By explicitly calibrating alternative metrics across a wide density gradient throughout most of the species' distribution, this study provides a transferable methodological framework for estimating density in low-density wildlife populations and the first continent-wide, standardized density assessment of a carnivore species. From a management perspective, our findings identify populations that may be most vulnerable, particularly those with the lowest densities, and highlight the need to prioritize absolute abundance monitoring.

An outbreak of New World screwworm has recently been spreading across Mexico, after more than 30 years of absence. The sterile insect technique, which consists of the massive release of sterilized males, has proven to be one of the most efficient methods for controlling the screwworm pest. However, given the limited number of sterile males available, improving the release strategy is critical. We propose a mathematical model of population dynamics adapted to the biology of Cochliomyia hominivorax and derive a feedback control function to determine the number of sterile males to release. We further construct a Luenberger observer to estimate wild fly populations from infected animal counts--the variable monitored by Mexican sanitary authorities--enabling field implementation of the control function. We show that eradication is achievable within approximately 60-100 weeks and that eradication time is governed primarily by the intrinsic biology of the system rather than by infestation magnitude. We then extend the model to a spatially explicit framework and show that when sterile male releases are applied at the outbreak focus and within a 120 km radius, eradication of the pest is attainable.

ContextThe movement of vertebrate hosts across urbanized landscapes can play a key role in the transmission of direct-contact diseases. Understanding how wildlife hosts move in urban landscapes, and how transmission is affected by their landscape-constrained and disease-altered movements, is imperative for better predicting the spread of disease. ObjectiveWe assess how the movement of red foxes (Vulpes vulpes) according to landcover type, and their infection status, affect the spread of mange (caused by Sarcoptes scabiei) in an urbanized landscape. MethodsWe developed a mange transmission model (MTM) using an agent-based model to compare two movement behaviours of foxes in Scarborough (Ontario, Canada): random and landcover-based. We further assessed the effects of movement on disease transmission by considering the foxs infection status and comparing a range of movement probability scenarios. We quantified the number of effective contact events and the effective reproduction number (Re) according to each scenario. ResultsWe found that both landcover-dependent movement and infection status influenced the spread of mange within fox populations. The number of effective contact events and effective reproduction number Re was greatest when landscape heterogeneity was included in the model and foxes moved through paths of least resistance to movement, and when susceptible and infected foxes had an equal probability of leaving a fragmented habitat patch. ConclusionsOur findings suggest that mange spread may be accelerated along movement corridors in fragmented, heterogenous landscapes. As urban areas expand and remnant habitat within these is further lost and animals are relegated to fewer movement pathways, disease transmission may increase.

The biodiversity crisis is exacerbated by persistent gaps in taxonomic, geographic, and conservation knowledge. This study provides a comprehensive assessment of Linnean, Wallacean, and conservation shortfalls in Rivulidae (Cyprinodontiformes), one of the most diverse families of Neotropical freshwater fishes. To this end, an extensive dataset was compiled, comprising 494 valid species, 51 synonyms, 3,419 occurrence records, and information on life cycle, distribution, and conservation status. The Linnean shortfall remains open: after 1975, the rate of species description increased sharply, and estimates indicate that 103 species remain undescribed (95% CI: 53-212). The year of species description was influenced by detectability and accessibility factors, with larger species, more widely distributed species, and species occurring in more densely populated areas being described earlier. The Wallacean shortfall was broad and spatially uneven: only 5.34% of the region was considered adequately sampled. The conservation shortfall was also substantial: 179 species are threatened with extinction, 69 species remain Not Evaluated, and 62 are Data Deficient. The mean time between taxonomic description and first IUCN extinction-risk assessment was 25.7 years. Moreover, 59.9% of species have no records within protected areas, including 69.3% of threatened species. These findings synthesize an urgent challenge: biodiversity knowledge and conservation shortfalls must be overcome simultaneously to protect species that are still being discovered, remain poorly documented spatially, and are restricted to habitats under intense anthropogenic pressure. The conservation of Rivulidae cannot wait for complete knowledge; action amid uncertainty is necessary to prevent both known and unknown species from disappearing.

Despite the ecological importance of ungulate migrations, we lack a complete understanding of why some ungulates migrate and others do not. Though progress has been made towards understanding differences across species and between populations, migratory behavior varies even within populations: in many populations, some individuals remain behind as residents (partial migration). Theoretical population-level work has suggested that these different migratory tactics can coexist, but such approaches stop short of providing insights into how individuals make the decision to stay or go each year. Using long-term data from three ungulate populations, we find that individuals probabilities of migrating are highly variable across years, which points to a non-trivial context-dependent decision-making process, whose underlying mechanisms must be probed via individual-level modeling. Drawing on existing knowledge, we propose a decision-making model of ungulate migration onset wherein individuals probabilistically decide to start migrating based on the local intensity of environmental and/or social cues. Residents arise as a robust collective organization phenomenon in our model. At sufficiently large population sizes, the number of residents is invariant with total population size, consistent with empirical patterns. Instead, resident numbers are influenced by the severity of the bad season, by relevant character differences among individuals, and by how individuals contribute and respond to environmental and/or social cues; for instance, when social cues contribute to decision-making in addition to environmental ones, fewer residents result, and migration is more likely to be complete. Overall, our model provides a potential mechanistic explanation for how residents might emerge within migratory ungulate populations.